Insight into the Dynamics of Fractional Maxwell Nano-Fluids Subject to Entropy Generation, Lorentz Force and Heat Source via Finite Difference Scheme

In recent times, the loss of useful energy and solutions to those energy challenges have a wide scope in different areas of engineering. This work focuses on entropy analysis for unsteady viscoelastic fluids. The momentum boundary layer and thermal boundary layer are described under the effects of a magnetic field in the absence of an induced magnetic field. The study of a fractional model of Maxwell nanofluid by partial differential equation using Caputo time differential operator can well address the memory effect. Using transformations, the fractional ordered partial differential equations (PDEs) are transfigured into dimensionless PDEs. Numerical results for fractional Maxwell nanofluids flow and heat transfer are driven graphically. The Bejan number is obtained following the suggested transformation of dimensionless quantities like entropy generation. A mathematical model of entropy generation, Bejan number, Nusselt number and skin friction are developed for nanofluids. Effects of different physical parameters like Brickman number, Prandtl number, Grashof number and Hartmann number are illustrated graphically by MAPLE. Results depict that the addition of nanoparticles in base-fluid controls the entropy generation that enhances the thermal conductivity and application of magnetic field has strong effects on the heat transfer of fractional Maxwell fluids. An increasing behavior in entropy generation is noticed in the presence of source term and thermal radiation parameter.


Introduction
The best way to represent the natural phenomena is by using differential equations (DEs) with suitable boundary conditions. Recently, the fractional (non-integral) order of DEs has gained much interest because of their vast scope in many engineering fields like food engineering, oceanography, chemical reactions and chaos. All fractional derivatives like Riemann Levillie, Caputo, Caputo Fabrizio, Antangna Beleanu are used widely in practice [1][2][3], as a deep study can easily be handled using this approach. Natural convection along a vertical wall and cylinder has been explained using Caputo time-fractional derivatives and fractional derivatives [4,5]. The significant difference between fractional and ordinary fluid flow emerges at different times. Furthermore, with large time values, increasing the fractional parameter's value increases velocity. Fractional calculus covers the complex structure of viscoelastic fluids in various research areas of glass fiber production, exotic lubricants, colloidal solutions, extraction of polymer solutions and cooling processes [6]. A simplified Phan-Thien-Tanner model for viscoelastic fluids was investigated 2 of 18 analytically in [7] with different physical parameters. It was established by Sheikh, N.A. et al. [8] that it is not enough to get experimental data by using a conventional derivative model for Maxwell fluids rather than fractional operators. Moreover, it was discovered that a link exists between Maxwell's constitutive equation and molecular theory [9]. Few other investigations for analytical results have been done in many research articles. But in the rheological perspective, viscoelastic fractional order models have attracted much interest due to their wide range of applications. Fractional Maxwell models are derived by substituting the conventional derivative in the known Maxwell model of stress-strain expression by fractional-order derivatives. Analytical solutions of fractional Maxwell, fractional generalized Maxwell model, fractional second grade, third grade, Oldroyd-B model etc., via analytical techniques of Laplace transform, Fourier transform, Weber transform and Hankel transform are obtained in [10][11][12][13][14].
Many researchers focus on such physical processes involving entropy generation-EG. It is a well-known fact that all physical problems, especially heat transfer, involve entropy generation-EG. Entropy generation-EG plays a vital role in fluid dynamics. The first law of thermodynamics moves around heat transfer processes, whereas the second law of thermodynamics moves around the entropy generation of the system. Entropy generation-EG tells the feasibility and efficiency of the system. In other words, entropy describes the ways a process can control energy loss.
Since entropy generation-EG is important and happens in almost all thermo-dynamical processes, many researchers have worked in this direction. Among the various researchers who have done splendid work on entropy generation-EG, Adrian Bejan, has published many articles and books [15][16][17][18][19][20][21]. Considerable work is done by Pranab Kumar Mondal [22] for irreversibility analysis of Couette flow while applying weak and relatively strong pressure gradient. The fluid dynamics are studied by variation of volumetric entropy generation number and Bejan number.
To control the entropy generation-EG, addition of nanoparticles in base fluid gave a new direction to heat transfer problems (the formation of nanofluids). Nanofluids are nano-sized particles accomplished with base-fluids, i.e., Common nanoparticles Cu, Ag, Au and Fe. Water, engine oil and ethylene glycol are common base fluids. To enhance the thermal conductivity of base fluids, the idea of nanofluids is given by Choi et al. in [23] for the first time. Later on, Tiwari and Das discussed the effects of different shapes of nano-sized particles on thermal enhancement [24]. Using square shape cavity, different aspects of ballshaped, cylindrical-shaped, and rod-shaped nanoparticles were studied [25]. Analytical results for the effects of nanofluids were driven by [26]. In [27], analytical results are investigated for temperature profile and Nusselt number under the effects of viscous dissipation and porous media. Using a traditional approach, numerical results were obtained, expressing the radial and tangential momentum across the disk decreases for higher Lorentz forces and slip factor [28].
Like many other key sources, entropy generation-EG includes viscous dissipation, chemical reaction, heat and mass transfer, heat convection and conduction and electrical conduction. Considerable work by B. Mahathish et al. has been done via the spectral quasilinearization (SQL) method for entropy analysis and solving a Williamson model, which provides a base for fractional models. An analytical approach is applied to investigate the effects of viscous dissipation and limiting effects of Nusselt for temperature profile [29]. Bejan number-BN is the ratio of entropy generation-EG due to heat and total entropy generation-EG of system. But Awed in [30] gave a new definition to Bejan number-BN. Bejan number-BN describes the effects of magnetic field irreversibility and fluid friction irreversibility. An investigation has been conducted on laminar falling liquid film along an inclined heated surface in [31]. In [32][33][34] and many other articles, the exact and numerical solution of entropy generation-EG is published. B. Mahathish numerically investigates the effects of quadratic variation of density-temperature (quadratic convection) and the quadratic Rosseland thermal radiation using the modified Bongiorno Model (MBM) [35]. Similar analytical as well as numerical approaches can be seen in [36][37][38][39] for such different viscoelastic models.
Numerical investigations and many others were used in [40][41][42][43]; however, various research gaps are found in these articles, which are still not addressed. Such as how heat transfer can be enhanced by adding the nanoparticle to the base fluid? What is the effect of Lorentz forces on flow dynamics? What do numerical results predict about Skin friction and Nusselt number? What are the formulation of entropy generation and Bejan number in the presence of thermo-physical properties of nanoparticles? In reviewing these research gaps in literature surveys, the main task of this article is to develop a mathematical model of momentum and heat for fractional Maxwell nano-fluids for detailed insights into Lorentz force, heat source/sink, and Nusselt number and entropy generation-EG. A new definition of Bejan number-BN is introduced by adding the coefficients of thermo-physical properties of nanoparticles. Water is taken as base fluid, whereas Cu and Al 2 O 3 are the nanoparticles used for graphical results of velocity and temperature profiles. The problem is first modeled fractionally by applying the definition of Caputo time derivative, then using transformation; a dimensionless analysis is done.

Mathematical Model and Formulations
Considering magnetohydrodynamic (MHD) flow of incompressible and unsteady fluid along the infinite vertical plate. The induced magnetic field and pressure gradient are neglected. Initially, at time t = 0, fluid has velocity zero and has a constant temperature θ ∞ . With the passage of time temperature of the system rise to θ w . The fluid flow is considered along x − direction. The magnetic field is applied in the y − direction as illustrated in Figure 1 below. (MBM) [35]. Similar analytical as well as numerical approaches can be seen in [36][37][38][39] for such different viscoelastic models. Numerical investigations and many others were used in [40][41][42][43]; however, various research gaps are found in these articles, which are still not addressed. Such as how heat transfer can be enhanced by adding the nanoparticle to the base fluid? What is the effect of Lorentz forces on flow dynamics? What do numerical results predict about Skin friction and Nusselt number? What are the formulation of entropy generation and Bejan number in the presence of thermo-physical properties of nanoparticles? In reviewing these research gaps in literature surveys, the main task of this article is to develop a mathematical model of momentum and heat for fractional Maxwell nano-fluids for detailed insights into Lorentz force, heat source/sink, and Nusselt number and entropy generation-EG. A new definition of Bejan number-BN is introduced by adding the coefficients of thermo-physical properties of nanoparticles. Water is taken as base fluid, whereas and are the nanoparticles used for graphical results of velocity and temperature profiles. The problem is first modeled fractionally by applying the definition of Caputo time derivative, then using transformation; a dimensionless analysis is done. The resulted model is solved by a numerical technique of finite difference scheme. Plots are drawn for Bejan number-BN against the Brickman number, the square of Hartmann number and Ω dimensionless temperature difference. Moreover, some graphical results are extended to evaluate the Nusselt number and Skin friction . These results are computed via mathematical software MAPLE.

Mathematical Model and Formulations
Considering magnetohydrodynamic (MHD) flow of incompressible and unsteady fluid along the infinite vertical plate. The induced magnetic field and pressure gradient are neglected. Initially, at time = 0, fluid has velocity zero and has a constant temperature . With the passage of time temperature of the system rise to . The fluid flow is considered along − direction. The magnetic field is applied in the − direction as illustrated in Figure 1 below; Taking into account the Boussinesq approximation, the assumptions of the system are, • Flow is unsteady, incompressible and 1-dimensional.
• Pressure gradient is absent.
• Body force is significant.
• Magnetic field is applied (ignoring the induced magnetic field). Taking into account the Boussinesq approximation, the assumptions of the system are,

•
Pressure gradient is absent. • Body force is significant. • Magnetic field is applied (ignoring the induced magnetic field).
Then the equation of continuity is restricted and takes the following form (can be seen in [44]); ∇·V = 0 (1) But the Navier-Stoke equation [45,46] takes the form, where ρ n f , T, J, B, g, β n f , θ and θ ∞ are dynamic viscosity of nanofluid, Cauchy stress tensor, current density, total magnetic field, gravitational acceleration, thermal expansion coefficient, the temperature of nanofluid and ambient temperature, respectively. The stress tensor for Maxwell fluids in [47] as, where In these expressions S, A 1 , λ 1 , µ, D/Dt represents the extra stress tensor, Rivline-Ericksen tensor, kinematic viscosity and material time derivative. Also A 1 is expressed in [48] as; And δS δt By following the assumption of the problem, the Maxwell equation can be written in [49] as; Keep in view that B = B 0 + b 0 is the sum of applied and induced magnetic field (neglected). Using Equations (3)- (7), the Equation (2) takes the form; Multiplying both sides of Equation (8) But the fractional constitutive equation for Maxwell fluids is given by [50], This constitutive relation contains D α t , Caputo fractional operator is defined in [51] as; With Γ(·) as the Gamma function defined in [51] by; Using Equation (10) into Equation (9) 1 Nanomaterials 2022, 12, 1745 5 of 18 The first law of thermodynamics [46,52] is; In this equation ρC p n f , K n f , q r and Q are constant of heat capacity, coefficient of thermal conductivity, radiative heat flux, and nanofluid thermal conductivity, respectively.
By using the Rosselands approximation for fluids that are considered optically thick, the radiative heat flux q r expressed in [53] is given as; In this expression σ * and k * are the Stefan-Boltzmann constant and mean spectral absorption constant, respectively. Approximating θ 4 by a Taylor's series expansion in the neighborhood of θ ∞ and neglecting higher power.
∞, (neglecting higher power), so radiative heat flux is Then Equation (12) becomes The second law of thermodynamics is given by [46], In which E G and σ n f are volumetric local entropy generation and electrical conductivity. The proposed boundary and initial conditions of this physical phenomenon are defined in [54] below; Using the following transformation and thermophysical properties of nanoparticles (can be seen in) [55], the dimensionless governing equation for velocity and temperature profile are obtained; Nanomaterials 2022, 12, 1745 6 of 18 The dimensionless velocity and temperature profile of the problem is given, and after omitting ( * ) notation for the sack of brevity of mathematical modeling.
Moreover, the temperature equation takes the following form In Equations (22) and (23) b 1 , b 2 , b 3 and b 4 are the ratio of thermophysical properties given by Additionally, Gr, Ha, Pr, Nr and Q 0 are the Grashof number, Hartmann number, Prandtl number, radiation parameter and heat generation parameter, respectively, defined in [45] as; Additionally, the non-dimensional initial and boundary conditions are The non-dimensional governing equation for velocity and temperature profile in Equations (22) and (23), with dimensionless initial and boundary conditions in Equations (24) and (25), represents the unsteady, incompressible flow fractional Maxwell nanofluids phenomena under the influence of magnetic fields. Water is taken as base fluid, but Cu or Al 2 O 3 are the nanoparticles considered for nanofluid preparation. For the numerical results the following Table 1 containing thermo-physical properties of nanoparticles and base fluid will be utilized.

Numerical Procedure
The finite difference scheme is a very efficient and powerful tool to investigate the numerical solutions of the problem arising and mathematical physics and mechanics. In this context, this section is dedicated to extending the finite-difference scheme to tackle the obtained set of fractional-order fluid models and heat transfer. For this, the discretization of the derivative of fractional-order of u, u t and u yy are specified as, Nanomaterials 2022, 12, 1745 and the nonlinear term is approximated by means of the following concept where, d α l = −l 1−α + (1 + l) 1−α for l = 1, 2, 3, . . . , j. Now, the rectilinear grid is assumed to examine the solution of the governing set of fractional-order fluid problems and heat transfer having grid spacing ∆y > 0, ∆t > 0 in the direction of space and time separately, where ∆t = T N , ∆y = L M for ∆y, ∆t from Z + . The inner points y i , t j in the discussed domain Ω = [0, T] × [0, L] are given as i∆y = y i and j∆t = t j . The discretization of the governing set of fractional-order fluid problems and heat transfer at y i , t j is given as, for j = 1, 2, 3, . . . , N − 1, i = 1, 2, 3, . . . , N − 1, with the following initial and boundary conditions,

Numerical Analysis and Discussion
Test Problem. For the validation of the applied scheme, a test problem is considered as In the discussed problem, the conditions are given below, and the source term can be selected against the choice of fractional-order derivative.
u(y, 0) = u t (y, t) = u(∞, t) = 0 and u(0, t ) = e at Since such a physical problem contains u(y, t) = y(y − t)t 2 as the exact solution. Its accuracy has been checked by a number of simulations for the proposed scheme. The plots in Figure 2a,b are drawn for the maximum absolute error (MAE) and computational order of convergence (COC) for different ranges of N, which is N = 10, 20, 40, 80, 160, 320, 640.
selected against the choice of fractional-order derivative.
The convergence of the applied scheme is observed against the selection of each fractional-order derivative, and its convergence order enhances for → 1. Figure 2c,d contains the L∞-norm between consecutive solutions that is − and | − | when 0 ≤ i, j ≤ N, M = 500. Again, it is found that the proposed scheme is very efficient, accurate and reliable for this problem. It also demonstrates that the solution is stable against the selection of fractional order and mesh parameters. The convergence of the applied scheme is observed against the selection of each fractional-order derivative, and its convergence order enhances for α → 1 . Figure 2c,d contains the L ∞ -norm between consecutive solutions that is u j+1 − u j ∞ and |u i+1 − u i | ∞ when 0 ≤ i, j ≤ N, M = 500. Again, it is found that the proposed scheme is very efficient, accurate and reliable for this problem. It also demonstrates that the solution is stable against the selection of fractional order and mesh parameters.

Entropy Generation
For viscous fluid flow in a magnetic field, the volumetric rate of local entropy generation E G is defined in [57] as; That is the sum of entropy generation due to heat transfer, due to fluids friction and due to the magnetic field effect, separately mentioned here (can be seen in [33]) Combining all the results The dimensionless entropy generation calculated is; where E 0 , is the characteristic entropy generation rate, Br is the Brickman number, H 2 a = M is the square of the Hartmann Number and Ω is the dimensionless temperature difference.

Bejan Number
The Bejan number is the irreversibility distribution parameter which is expressed mathematically as (can be seen in [58]); Be = Entropy generation due heat transfer Total Entropy generation , That is Br Ω ∂u ∂y 2 + a 6 Br Ω H 2 a (u) 2 Then the reduced expression for the Bejan number is given as; Br Ω ∂u ∂y Br With b 5 = a 3 a 5 , b 6 = a 6 a 5 .

Skin Friction and Nusselt Number
For measuring shear stress and heat transfer effects in an ordinary integer order system, local skin friction and Nusselt number are defined in [59] as; and The skin friction coefficient and local Nusselt number for (FMF) can be written by using Equation (5), that is, the fractional stress tensor for Maxwell fluid on the plate with fractional time Caputo derivative (detail can be seen in [60]).
The non-dimensional form of Equations (38) and (39) is given as

Results and Discussion
This section of the article deliberates the detailed results and discusses the plots driven against different physical parameters like H a , Pr, Nr, Gr, φ, α, β and Q 0 representing the magnetic field parameter (the square of Hartmann Number), Prandtl number, radiation parameter, Grashof Number, volumetric fraction of nanoparticle, fractional order parameters and heat generation parameter, respectively. The behavior of these aforementioned physical parameters on dimensionless velocity profile u(y, t), temperature profile θ(y, t), Bejan number-BN Be, Skin friction S f and Nusselt number Nu are drawn by mathematical software MAPLE. The mathematical fractional model of Maxwell nanofluid is developed by using the Caputo time fractional operator. After transforming the governing equations to a dimensionless governing model, the finite difference method (FDM) is used for the discretization of the model. FDM is a strong tool for dealing with such kinds of problems. The simulation is performed by developing and executing codes.
Results are obtained by solving Equations (22), (23), (32), (39) and (40) with initial and boundary conditions illustrated in Equations (24) and (25) and physical properties of nanoparticles in Equation (21) and Table 1. Various suitable ranges of physical paramet (H a = 1, 2, 5), (Pr = 3.5, 6.2, 15), (Gr = 0, 1, 2), (φ = 0.01, 0.1, 0.2), (α = 0.2, 0.4, 0.6, 0.8, 1), (Nr = 0, 2, 5), (Q 0 = 0, 2, 5) for heat transfer, velocity analysis, skin friction, Nusselt number, entropy generation and Bejan number-BN are considered, and also particular exertion has been done on the effects of these parameters for heat enhancement. Figure 3 shows the effect of the magnetic field parameter on the velocity profile of fractional Maxwell fluids. It can be seen that velocity decreased by increasing the value of H a because of the Lorentz's force. The rise in H a caused strengthening in the Lorentz force, which increases the internal resistance to flow particles; consequently, fluid velocity decreased. Whereas an increase in Skin friction occurs, as shown in Figure 3b. Since Bejan number Be is the ratio of total entropy generation to entropy generation due to heat transfer. Thus, advancement in H a boosts the Bejan number-BN Be and can be noticed in Figure 3c. An opposite behavior of total entropy generation N S relative to Bejan number-BN Be appears against the increasing value of H a , as can be seen in Figure 3d. Figure 4 depicts the effects of the Grashof number Gr on velocity profile u(y, t). The enhancement in value of Gr results in increasing the fluid velocity, which can be physically justified as the increasing value of Gr means lower the viscous forces and hence increasing the velocity of fractional Maxwell fluids. While Skin friction S f increases, the increasing the value of Gr can be noticed in Figure 4b. A slight effect of Gr on Bejan number-BN Be and entropy generation N s can be seen in Figure 4c,d, respectively. Similar results are reported by Sarojamma, G., et al. in [61] for comparison. Gr is the ratio of buoyancy forces to viscous force, increasing value of Gr results in laminar boundary layer and vice versa. Because high value of Gr give rise to the temperature of molecules, consequently, the intermolecular forces become weak. Thus, velocity profile u(y, t) also rises. On the other hand, fluid particles collectively gain momentum as Gr increases, so additional heat is lost nearby, that is why skin friction reduces as shown in Figure 4b. Since Bejan number Be is ratio of entropy generation due to heat to the total entropy generation of the system, that is why Gr reduces the value of Bejan number-BN Be deliberated in Figure 4c.  increases, the increasing the value of can be noticed in Figure 4b. A slight effect of on Bejan number-BN and entropy generation can be seen in Figures 4c and 4d, respectively. Similar results are reported by Sarojamma, G., et al. in [61] for comparison.
is the ratio of buoyancy forces to viscous force, increasing value of results in laminar boundary layer and vice versa. Because high value of give rise to the temperature of molecules, consequently, the intermolecular forces become weak. Thus, velocity profile ( , ) also rises. On the other hand, fluid particles collectively gain momentum as increases, so additional heat is lost nearby, that is why skin friction reduces as shown in Figure 4b. Since Bejan number is ratio of entropy generation due to heat to the total entropy generation of the system, that is why reduces the value of Bejan number-BN deliberated in 4c.  Figure 5 shows the effects of the volume fraction parameter φ on fluid velocity, and it is examined that advancement in controlled volume fraction parameter φ lowers the velocity profile due to the effect; that is, the viscosity of fractional Maxwell fluids increases by increasing φ. On the other hand, the addition of nanoparticles in the base fluid causes an improvement in the heat transfer rate at the boundary layer. As the thermal conductivity of base fluid is enhanced, consequently increasing the fluid's internal temperature can be seen in Figure 5b. Since Nusselt number Nu is the ratio of convective heat transfer to conductive heat transfer, increasing φ has a decreasing relation with Nusselt number Nu. This is because the skin friction decreases with the passage of time against the volume fraction parameter φ. These results for skin friction and Nusselt number Nu are viewed in Figure 5c,d, respectively. An increasing trend is noticed for the Bejan number Be and can be seen in Figure 5e as the heat transfer rate becoming better by increasing φ. Whereas in Figure 4f Figure 5 shows the effects of the volume fraction parameter on fluid velocity, and it is examined that advancement in controlled volume fraction parameter lowers the velocity profile due to the effect; that is, the viscosity of fractional Maxwell fluids increases by increasing . On the other hand, the addition of nanoparticles in the base fluid causes an improvement in the heat transfer rate at the boundary layer. As the thermal conductivity of base fluid is enhanced, consequently increasing the fluid's internal temperature can be seen in Figure 5b. Since Nusselt number is the ratio of convective heat transfer to conductive heat transfer, increasing has a decreasing relation with Nusselt number . This is because the skin friction decreases with the passage of time against the volume fraction parameter . These results for skin friction and Nusselt number are viewed in Figures 5c and 5d, respectively. An increasing trend is noticed for the Bejan number and can be seen in Figure 5e as the heat transfer rate becoming better by increasing . Whereas in 4f entropy generation enhanced by increasing the value of . Because the viscosity of fluid increased by increasing . In these following plots in Figure 6, results are drawn for velocity profile u(y, t), temperature profile θ(y, t), skin friction S f , Nusselt number Nu and Bejan number-BN Be against fractional parameter α.
The reason is that the gradual increase in fractional parameter α gives rise to the viscosity of the nanofluid. This means intermolecular forces between the nanoparticles and base-fluid particles increase; consequently, Brownian motion of particles reduces, that is, a decrease in velocity profile u(y, t) occurs. The consequences of fractional order parameter α for velocity profile u(y, t) are inverse as depicted in Figure 6a. On the other hand, α have direct relation for θ(y, t), as shown in Figure 6b. Figure 6c shows that for an increasing value of α the entropy generation N s decreases, whereas a rise in α varies Be directly, deliberated in Figure 6d. Figure 7 shows the effects of thermal radiation parameter N r on temperature θ(y, t) of the fluid. Applying thermal radiations gives rise to the temperature of particles of nanofluids. Hence, the particles' kinetic energy increases, and the rate of collision between the particles of the nanofluids becomes high, which is why a rising increase in temperature profile occurs. Therefore it is concluded that increasing the value of N r causes an increase in fluid temperature. Nu is the ratio of convective to conductive heat transfer across the boundary, but enhancement in conduction occurs with the addition of nanoparticle therefore, a decrease occurs in Nu, depicted in Figure 7a,b, respectively. A decrease in Bejan number-BN Be can be noticed in Figure 7c, whereas entropy decreases initially and then increases gradually, as can be seen in Figure 7d. Similar results are reported in [61].
Prandtl number Pr is the dimensionless number and is the ratio of momentum to thermal diffusivity. It is a fluid property but does not have any dependence on flow type. Thus, an increase in Pr means heat transfer is favored to occur by momentum, not conduction. This parameter controls the relative thickness and thermal boundary layer in heat transfer problems. Lowering the value of Pr means the heat diffuses spontaneously as compared to momentum, which thickens the thermal boundary layer rather than the momentum boundary layer. Therefore, an increase in Pr decreases the temperature profile θ(y, t) of fractional Maxwell fluids as expressed in Figure 8a, which expectedly decreases the Nusselt number Nu. Since Bejan number Be has an inverse relation with entropy generation N s , due to heat transfer, that is, Be decreases with the increase in value of Pr as illustrated in Figure 8. Figure 9 shows the effects of the heat source term Q 0 on temperature profile θ(y, t), since heat source gives rise to temperature profile as shown in Figure 9a, but with the passage of time, it is noticed that the Nusselt number Nu decrease is dependent inversely on the conduction of heat. Thus, entropy generation N s expectedly decreases as shown in Figure 9d. This phenomenon gives rise to the Bejan number Be because entropy generation due to heat transfer has an inverse relation with the Bejan number Be, as shown in Figure 9c. In these following plots in Figure 6, results are drawn for velocity profile ( , ), temperature profile ( , ) , skin friction , Nusselt number and Bejan number-BN against fractional parameter .
The reason is that the gradual increase in fractional parameter gives rise to the viscosity of the nanofluid. This means intermolecular forces between the nanoparticles and  Figure 7 shows the effects of thermal radiation parameter on temperature ( , ) of the fluid. Applying thermal radiations gives rise to the temperature of particles of nanofluids. Hence, the particles' kinetic energy increases, and the rate of collision between the particles of the nanofluids becomes high, which is why a rising increase in temperature profile occurs. Therefore it is concluded that increasing the value of causes an increase in fluid temperature.
is the ratio of convective to conductive heat transfer across the boundary, but enhancement in conduction occurs with the addition of nanoparticle therefore, a decrease occurs in , depicted in 7a and 7b, respectively. A decrease in Bejan number-BN can be noticed in Figure 7c, whereas entropy decreases initially and then increases gradually, as can be seen in Figure 7d. Similar results are reported in [61]. Prandtl number is the dimensionless number and is the ratio of momentum to thermal diffusivity. It is a fluid property but does not have any dependence on flow type. Thus, an increase in means heat transfer is favored to occur by momentum, not conduction. This parameter controls the relative thickness and thermal boundary layer in heat transfer problems. Lowering the value of means the heat diffuses spontaneously as compared to momentum, which thickens the thermal boundary layer rather than the mo-  Figure 9 shows the effects of the heat source term on temperature profile ( , ), since heat source gives rise to temperature profile as shown in Figure 9a, but with the passage of time, it is noticed that the Nusselt number decrease is dependent inversely on the conduction of heat. Thus, entropy generation expectedly decreases as shown in Figure 9d. This phenomenon gives rise to the Bejan number because entropy generation due to heat transfer has an inverse relation with the Bejan number , as shown in Figure 9c.

Conclusions
The graphical analysis of fractional Maxwell nanofluids is made in this article under the influence of a magnetic field (ignoring induced magnetic field). The pressure gradient is supposed to be absent. Effects of different physical parameters are drawn by using the mathematical software MAPLE. The model is formulated by applying the Caputo time

Conclusions
The graphical analysis of fractional Maxwell nanofluids is made in this article under the influence of a magnetic field (ignoring induced magnetic field). The pressure gradient is supposed to be absent. Effects of different physical parameters are drawn by using the mathematical software MAPLE. The model is formulated by applying the Caputo time derivative. Using suitable transformations, governing equations are made dimensionless.
MHD fractional Maxwell nanofluids are studied numerically; quantities like entropy generation, Bejan number, Skin friction, and Nusselt number are investigated using the finite difference method.
Hence key findings of this study are given below; a-It is noted that for increasing the value of fractional order parameter α, the velocity profile u(y, t) decreases, whereas the temperature profile θ(y, t) increases. b-The addition of nanoparticles to base fluid enhances the thermal conductivity of fractional Maxwell nanofluids, increasing the value of volume fraction of nanoparticles φ and decreasing entropy generation N s . c-The magnetic field effect influences the temperature θ(y, t) and velocity u(y, t) profile with inverse and direct behavior, respectively. d-Nusselt number increases with the variation in Pr, and a decrease occurs in Nu with the increase in thermal radiation parameter. e-The temperature profile varies directly with the thermal radiation parameter Nr, and increasing the value of Nr decreases the Nusselt number. Whereas entropy generation N s increases, and the Bejan number decrease with a rising value of Nr.
The solution obtained via the finite difference method is excellent in agreement with the test problem and existing results which shows that the finite difference method is a strong and reliable technique to deal with such kind of complex models and it gives a key direction for further study.